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Abstract 

Using the statistical inference method, a non-relativistic, spinless, non-hnear 
quantum dynamical equation is derived with the Fisher information metric substi- 
tuted by the Jensen-Shannon distance information. Among all possible implications, 
it is shown that the non-linear Schrodinger equation preserves the symplectic struc- 
ture of the complex Hilbert space, hence a Hamiltonian dynamics. The canonically 
projected dynamics is obtained on the corresponding projective Hilbert space of 
pure state density operators. 



^ j 1 Introduction 



The time-evolution equation of quantum states, namely the linear Schrodinger equation 
can be obtained by using the maximum uncertainty principle introduced by Jaynes [1], 
and using the statistical inference method [2, 3]. The idea, in essence, is to minimize the 
appropriate action functional, containing the statistical mechanical action and a term 
proportional to the Fisher information measure [4], and combine the two resulting real 
equations to one linear complex equation by introducing the wavefunction. 

Indeed, minimization of the Fisher information metric, to obtain the linear quantum 
theory, is subject to set of physical constraints. According to the Cramer-Rao inequality 
in statistical estimation theory [5], the variance (mean-square error) of any unbiased 
underlying parameter, 6, of a probability distribution, p{x;6), is maximized once the 
Fisher information measure is minimized, hence, the maximum entropy principle. 

The axiomatic approach in seeking generalizations beyond the linear theory is studied 
by Parwani [6]. In [7], he proposed a non-linear version of the Schrodinger equation 
based on the same information theoretic discussions by replacing the Fisher information 
measure with the relative entropy, or the Kullhack-Leihler information measure [8], in 
the appropriate action functional, which is partially motivated from the fact that for two 
close statistical distributions and to the lowest order, the latter measure is proportional 
to the former. 
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Applying the variational method would then lead to the new quantum Hamilton- 
Jacobi equation and the continuity equation, which again can be combined, by introduc- 
tion of the wavefunction, to form the non-linear Schrodinger equation. The non-linear 
term is a function of an invariant length scale which inspires the breaking of space-time 
symmetries, namely the rotational and Lorentz invariances in small scales, i.e. high 
energies. 

In the current paper, we aim at an optimized version of the non-relativistic, spinless, 
non- linear quantum dynamical equation, introduced in [7], which will be obtained by 
minimization of the bounded, smooth and always well-defined version of the relative en- 
tropy, called the Jensen-Shannon distance measure [9]. The common properties as well 
as the different features in comparison with the Kullback-Leibler non-linear quantum dy- 
namics will be studied. Moreover, the implications this rather information-geometrically 
motivated non-linear quantum dynamics may have on the structures of the linear theory 
and the dynamical properties on the complex Hilbert space will be investigated. 

The outline of the paper is as follows. In section 2, both Kullback-Leibler and Jensen- 
Shannon measures as well as the statistical inference method are reviewed and the metric 
property of the latter measure is highlighted in further details. In section 3, by using 
the statistical inference method, the non-linear quantum dynamics is obtained and the 
extension to many particles in higher dimensions and also its parametric versions are 
derived. In section 4, dynamical properties of the non-linear Schrodinger equation are 
studied, in the language of geometric quantum mechanics, followed by a brief discussion 
about the metric structure. Finally in section 5 we summarize the results. 



2 Information Measures and Statistical Inference Method 
2.1 Kullback-Leibler Information Measure 

The Kullback-Leibler (KL) divergence, or information measure, was first introduced in 
[8] as the amount of information for discrimination between two statistical distributions 
Pq{x) and pi{x). It is defined as 

Ikl{pi,Po) = [ pi{x)ln(^) d^xdt, (1) 

with po{x) as a reference probability distribution encoding a prior information about 
the ensemble. The probability distributions are the functions of collective continuous 
coordinates x = x^, with k = 1, 2, ... m. 

As it is evident from equation (1), KL information measure is asymmetric" with 
respect to two distributions, or in other words, I k l{po{x) , pi{x)) is not the same as 
I K l{pi{x) , Pq{x)) . Therefore, the KL measure cannot be considered as a distance along 
the curve connecting the points po ^-nd pi on the space of statistical distributions. 

Furthermore, the functional form of the equation (1) is not invariant under an arbi- 
trary re-par ametrization (except affine re-par ametrizations) of the integral (1), which is 
in contrast with the characteristics of the geodesies in Riemannian geometry being the 
invariant distance minimizing the curves under arbitrary re-parametrizations [10]. 
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Different symmetric versions of tlie KL divergence measure exist [11], according to 
various ways of taking the "average" of two KL measures, such as Jeffreys-Kullback- 
Leibler or J-divergence which was introduced in [12, 13] as 

J{po,Pi) = ^ [Ikl{po, Pi) + IkUpi, Po)] ■ (2) 

It can be considered as the average distance between two probabiUty distributions po and 
pi. The J-divergence is symmetric and also invariant under the re-parametrization but 
it doesn't satisfy the 'triangle' inequahty, and is therefore not considered as the metric 
distance in the space of statistical distributions. Other problems with the J-divergence 
are its 'unboundedness' and 'instability' in practical computations [14]. 
The following identity holds for the Ikl, 

Ikl{po,pi) + Ikl{pi,po)= / iF{e)de (3) 

Je=Q 

where p{x,9 = 0) = Pq{x) and p{x,6 = 1) = pi{x), and the Fisher information measure, 
Ip, is defined as [4] 

In other words, the integral of the Fisher information measure along the geodesic pg 
connecting po to pi is identified with the J-divergence. 

Choosing pi{x^) = p{x^ + Ax^), the KL measure and J-divergence can be expanded, 
and to their lowest orders one can obtain 

Ikl {pix' + Ax^),/,(x'=)) ^ -\ (Ax^y If{p{x)) (5) 

and 

j(p{x^),p{x^ + Ax'')) =-(Ax^)' If{p{x)) (6) 



over a Euclidean manifold {M,5ij) and 6^^ as the inverse of the Euclidean metric. 
2.2 Jensen-Shannon Information Measure 

To avoid the problems of Kullback-Leibler relative entropy, as stated in sub-section 2.1, 
the Jensen-Shannon information measure, Ijs, is introduced by Lin [9]. For a one- 
dimensional Euclidean parameter space it is defined as 

Ijs = ^{lKipo,Pi)+lK{pi,Po)) (7) 
where each Ik in the RHS of equation (7) is defined as 

Ik{po-,Pi)= I Po In ( ) dxdt (8) 

J \Po + PiJ 
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and it is called the K-divergence. 
Recalling the Shannon entropy 



H{p) = - p{x)ln{p{x))dx, (9) 



the Jensen-Shannon divergence can be written as 

Ijsipo, Pi) = H{ — ) — > 0. (10) 

As it is shown in [9], the following inequality illustrates the upper limit of the Jensen- 
Shannon information measure 

0<Ijs{po,Pi)<\jiPo,Pi)- (11) 

Similar to the discussion of the previous section, for two probability distribtions of p{x) 
and the shifted p{x + Ax), in the limit Ax — )• 0, one can find [15] 

Ijs^liAxflF. (12) 

A generalized version of the Jensen-Shannon divergence measure, for two probability 
densities with different weights, is also introduced in [9]. Assuming < tti, 7r2 < 1 as the 
weights for the probability densities pi and p2, with tti + tt2 = 1, the generalization of 
equation (10) is introduced as 

4s(/'l'P2) = H {TTipi+7r2P2) 

-TTlH ipi) - 7r2H {p2) > 0. (13) 

This is called the "vr-parametric Jensen-Shannon divergence". For vr = tti = 7r2 = 1/2, 
it reduces to the equation (7). In one dimension for simplicity, the vr-parametric K- 
divergence is then equal to 

iP (P1,P2) = /piln ( —) dxdt (14) 

and the Jensen-Shannon measure can be written as 

I^fkPl^P2) = \ {iP iPl,P2) + iP (P2,Pl)) . (15) 

In [7], the weight vr appeared as a regularising parameter, which in fact transforms 
the Kullback-Leibler information measure to the asymmetric K-divergence defined in 
equation (14), but not /j^- As in equation (12), in the limit Ax — t- and for two close 
probability densities p{x) and p{x + Ax) with the weights of vr and (1 — vr) respectively, 
we have 

{p{x),p{x + Ax)) (1 - vr) If (p(x)) . (16) 
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Therefore, according to [15, 16] and the equation (7), the Jensen-Shannon divergence 
is the symmetrized, well-defined and bounded version of the KL information measure 
satisfying the triangle inequality. Recalling the definition of a distant vector space, 
it turns out that the square root of the JS divergence denotes a metric in the space 
of statistical probability distributions [17, 18], i.e. the space {Ai, y/Ijs), with Ai as 
the set of statistical probability distributions, is a metric space. Using the results of 
Schoenberg [19], one can prove that the space of probability distributions, with the 
metric Djs{po, pi) = \/ Ijs{po, Pi), can be isometrically embedded in a real separable 
Hilbert space [20, 21]. The importance of the metric property of the Jensen-Shannon 
divergence will be discussed in the sub-section 2.3 below. 

2.3 Statistical Inference Method 

The dynamical equations of an ensemble of classical particles, moving in one dimension 
X for simplicity, can be obtained by minimizing the following classical action functional 

with respect to p and S as independent variables. The resulting equations are the classical 
Hamilton- Jacobi equation and the continuity equation. The latter equation implies that 
the total probability, J p{x) dx = c, is conserved. Therefore, in the space of statistical 
densities, the subspace of normalized densities, where c = 1 and p{x) > for all x € M, 
are selected during this process. 

To derive quantum dynamical equations, some physical constraints have to be im- 
posed to the above-mentioned procedure. As mentioned in section 1, according to the 
Cramer-Rao inequality, the variance (mean-square error) of any unbiased parameter, such 
as the particle's positions, is maximized if the Fisher information measure (4) is mini- 
mized. The constraint can be imposed by the method of Lagrange multipliers, namely 
by applying the variational technique to the modified action functional IcL + C^F, as 
in [2], and the resulting equations are the quantum Hamilton- Jacobi equation and the 
probability conservation equation. 

Geometrically, representing the space of statistical probability distributions as a man- 
ifold, the Fisher information measure defines the Riemmanian metric on it [22]. The 
embedding from the normalized submanifold of the statistical densities to the infinite- 
dimensional real Hilbert space of the square-integrable functions, L^(]R), is provided by 
the map p{x) — t- p{x) and the corresponding Riemannian metric on the embedded unit 
sphere of the real Hilbert space is the Fisher-Rao metric, which is proportional to the 
Fisher information metric in form of the log-likelihood functions [23, 24]. In the statisti- 
cal inference method to obtain the linear quantum dynamics [2], the polar representation 
of the wavefunctions, i.e. 'i/;(x) = ^ypexp{iS{x)/h), is used as the embedding map and 
one can show that the metric induced on the Lagrangian submanifold of the unit sphere 
in the complex Hilbert space, 7i, is in fact the Fisher-Rao metric, i.e. when dS = or 
S{x, t) = const [25]. 
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A rather physical axiomatic approach in seeking generahzations beyond hnear quan- 
tum theory, namely substituting the Fisher metric by the Kullback-Leibler information 
measure in the statistical inference method, is studied in [6]. Therefore, according to 
the axioms in [6] and the discussions of the sub-sections 2.1 and 2.2, the Jensen-Shannon 
metric (7) is an optimized measure in seeking generalizations beyond linear quantum the- 
ory. The construction of the non-linear quantum dynamics from other distance measures, 
such as Renyi's and Wooter's divergences, are also discussed in [7]. 

The significance of the distance or metric property of the Jensen-Shannon divergence 
is even more highlighted by recalling the notion of distinguishability between two pure 
quantum states in Wootters' paper [26], i.e. the angle between the associated rays in the 
Hilbert space is equivalent to the concept of statistical distance. More precisely, the met- 
ric property of the Jensen-Shannon divergence, i.e. Djs = V^JS^ of the statistical space 
M is vital, at least in the sense of preservation of the notion of distinguishability of pure 
states in the embedded submanifold of the Hilbert space, as a Kahler space. Moreover, 
the idea of distinguishability of pure states is required for the geometrical interpretation 
of the transition probabilities in the measurement process to remain meaningful, as the 
geodesic distance between two distinguished states in the (projective) Hilbert space, in 
the linear theory. 

In the subsequent section, by using the inference method of the maximum uncertainty 
principle, both non-parametric and parametric versions of the non-linear Schrodinger 
equation will be derived. 

3 Generalized Quantum Dynamics 
3.1 The Non-linear Schrodinger Equation 

The discussion of the previous sub-sections 2.2 and 2.3, motivates one to replace the 
Kullback-Leibler divergence by the Jensen-Shannon information measure introduced in 
equation (7), and minimize the action functional IcL -\- C^JS with respect to the real- 
valued functions p and S, as in [7]. 

Applying the variational method on Iql + C^JS^ with two shifted probability distri- 
butions pq = p and pi{x) = p{x + I) = pi and the Lagrange multiplier C, = h? /mP , the 
standard continuity equation and the modified quantum Hamilton- Jacobi equation are 
derived respectively as follows 




(18) 




(19) 



with the following quantum potential 



Qn = - In 



4p2 



(20) 



{P + Pi){P-i + P) 
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and p_i = p{x — I). Combining the two real equations through the mapping ijj = 
y/pexp{iS/h), the non-hnear complex Schrodinger equation is obtained as 



where the non-linear term is 



In 



{p + Pl)ip-l + p) 



The equations (20) and (22), are the altered versions of the quantum potential and 
the non-linearity term found in [7] respectively. Regarding the non-linear Schrodinger 
equation (21) and the modified quantum potential (20), the following remarks are in 
order: 

Remark 1. By expanding the Qf^f, one can easily verify that to the zeroth order of /, 
Qn — Q and so N{p) ~ 0, where Q stands for the ordinary quantum potential derived 
from the linear Schrodinger equation, i.e. the last term in the RHS of the equation (22). 
To the leading order, the expansion of the non-linear term N{p) is 

NlpM^-^C-^-'-^.O.Oi,) (23) 



64 m \ p2 p^ p* ^ 

where p' and p" denote dxP and d^p respectively. Therefore, the first leading order of 
expansion is of the order of t^. 

Remark 2. The scale invariance property of the non-linear term (22), namely N{'y'^p) = 
N(p), allows the solutions of the equation (21) to be normalizable, i.e. if tp is the solution, 
then so is the 7'i/'. If 7 is only a phase factor, then the solutions are globally (with respect 
to the configuration) phase invariant. 

Remark 3. For the Gaussian wavepackets, namely ^|) cx exp(— rr^/ft^), Qj\i leads to a 
repulsive quantum force for all x and it is the cause for wavefunction dispersion, as in 
the ordinary quantum potential. 

Remark 4. Since the ordinary continuity equation is preserved, one can imply that the 
quantum probability current j is also left unchanged. To show it explicitly, one can write 
the non-linear Schrodinger equation (21) and its complex conjugate for the generalized 

Hamiltonian H = Hq + N{p) and multiply them by ip* and —ip respectively. Using the 
fact that the non-linear term N(p) is a real- valued scalar functional, the summation of 
the obtained equations would lead to the standard continuity equation, with the quantum 
probability current j defined as in the linear quantum dynamics. It is then implied that 
the norm-squared is conserved, for pure quantum state ^ £ H of the non-linear 

Schrodinger equation (21). 
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3.2 Many Quantum Particles in Higher Dimensions 

The discussion of the previous section 3.1 can obviously be extended to higher dimensions 
with n particles evolving on a d-dimensional Euchdean physical space as the following. 
The classical action functional of equation (17) is modified, as in [7], as 

fdS , 5^^ dS dS , _\ 



for n particles, with x = x*, i,j = 1,. . . ,nd and the index (i) as the smallest integer 
> i/d, as in [7]. In other words, i = {k—l)d+l, . . . ,kd refers to the coordinates of the kth 
particle with the mass m^. Treating each generalized coordinate of the nd-dimensional 
configuration space separately, the extended Jensen-Shannon divergence (7) for the ith 
degree of freedom, i.e. Ijg{p,pi), is equal to 



P In — ' + In ' 



P + PiJ \P + P: 



dx' dt (25) 



where pi{x) = p{xi,X2, ■ ■ ■ , Xi + l[^, . . . , Xnd, t) and the symbol l^^ defined as [i] = i mod d 
with the convention /jg] = ^d, due to the fact that the invariant length scales l^^ are 
defined in the d-dimensional Euclidean physical space, but not in the nd-dimensional 
configuration space. 

As usual, one can apply the variational technique to the extended action functional 
IcL + C^JSi with the latter term defined as 

nd 

1=1 

and Q = combines the resulting continuity equation and the quantum Hamilton- 

Jacobi equation, through ^ = ^^exp{iS/h), to obtain the non-linear Schrodinger equa- 
tion as ^ 

ihdt^ix, t) = ( + V(x)] ^ + N(p)^ (27) 



with the non-linear term as 

nd 



8m(i) \ p p 



( _ ^^11^ ] (28) 



and p±i{x) = pixi,X2, ...,Xi± ^[j], . . .,Xnd,t). 

The non-linear Schrodinger equation (27) for n non-interacting particles decouples 
into n independent non-linear Schrodinger equations for each individual particle. In 
other words, when V {fi,r2, . . . , f„) = V^^\fi)i=BV^'^\f2) © • • •®V^'^\fn), the particle 
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is not influenced by an action on particle 'j' by change of V^^\fj), for either initially en- 
tangled or factorized (product) state of the composite system, i.e. ^ (ri, f2, . . . , r„; to) = 

Ip^'^^ in, to) ® (r2, to) ® • • • ^ V^"H^n, to). 

However, this is the direct consequence of the '^complete separability" of the non-linear 
term in equation (28), in the sense introduced in [27], which is defined as 

N{p) = Ni e N2 e • • • e A^„, (p^")) , (29) 

where p {fi,f2, ■ ■ ■ ,i^n) = P^^K^^i) ^ p''^H^2) ^ ■ ■ ■ ^ P^^\rn), with 



) = y • • • y 1^ (n, f2, . . . , fn;t)\^ f[ dfj, (30) 



as the reduced probability density of the ith particle, in which the integration for that 
particle is singled out. In this sense, the non-linear Schodinger equation for each particle 
is written independently as 

ihdti^^^ {U,t) = + iV.(/9«)) V'^'^ {ri,t) (31) 

where 

/?« = -^A, + y«(,-^). (32) 

As it is discussed in [27], a consistent extension from a single particle system to 
a many-particle composite system is guaranteed by the complete separability of the 
nonlinear Hamiltonian operator in terms of reduced densities of individual particles, i.e. 
the equation (29), which will not lead to the locality problem [28-30], for entangled or 
generally mixed state of the composite system. Furthermore, by using the approach 
introduced in [31], this construction can be generalized to the multiple-time correlation 
experiments on the many-particle composite system, without unphysical nonlocal effects, 
such as those proposed by Gisin in [32]. 

3.3 Parametric Non-linear Schrodinger Equation 

Recalling the definition of the vr-parametric Jensen-Shannon information measure in sub- 
section 2.2 and the discussion of the subsection 3.1, one is motivated to repeat the same 

(it) 

procedure, but this time with Ijg- For one particle moving in one dimension, adopting 
the weight vr for the the probability distribution p{x) and so (1 — vr) for the shifted density 
p{x + I) = pi, the action functional that has to be minimized is then Icl + V^^js- 

As in the previous sections, the Lagrange multiplier r] is fixed from comparing equa- 
tion (16) with /i^/8m, which is used in obtaining the linear quantum theory [2], and 
therefore 

^2 

''=47r(l-^)m/2- 
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The variational method with respect to p and S*, leads to the ordinary continuity equation 
and the parametrically modified quantum Hamilton-Jacobi equation with the parametric 
quantum potential term as the following 



(47r)-^n^ 
(1 — vr) mP 



In 



p(7rp+ (1 - tt)pi) 
(7rp_^ + (l-vr) 



(34) 



Combining the resulting parametric Hamilton-Jacobi and continuity equations, through 
the identification ^ = ^JJ)exp{iS/h), leads to the following 'parametric' non-linear 
Schrodinger equation, for one particle in one dimension for simplicity, as 



ih 



dip 
dt 



(35) 



with the non-linearity term defined as usual as 

N^^\p) = Q^;;^ -Q (36) 
where the Q is the ordinary quantum potential in the linear quantum dynamics. 



4 Dynamical Properties 
4.1 Background 

The geometrical settings of the linear quantum mechanics has been studied thoroughly so 
far, for instance in [33-38] and [39, 40] just to mention a few in the existing literature. It is 
well-known that the space of pure quantum states, namely the complex Hilbert space Ti, is 
equipped with the structure of a Kahler manifold, (Q, G, J) with as the non-degenerate 
symplectic form, G as the metric structure and J as the compatible complex structure, 
which is a linear operator representing multiplication by i and satisfying = —1. 

Considering "H as a real vector space, the Q and G can be identified with the real and 
imaginary parts of the Hermitian inner product [37, 39] 

{^\^) = ^G{^,<^)+'-n{^,^) (37) 
with the following relationship 

!^(^,$) = G(J^',$) (38) 

where (^|J$) = For both finite-dimensional and infinite-dimesional Hilbert 

space, one can show that the linear Schrodinger equation is equivalent to the Hamiltonian 
dynamics, generated by the expectation value of the associated Hamiltonian operator 
[33, 35, 37, 41], namely 

dHo{Y) = Q(^X^^,Y) (39) 
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for any vector field Y G XCH) and Hq = (^|-ffo^')/(^|^) and 

XfjJ^):=-^Ho^. (40) 

However, the Hilbert space Ti is not the true space of states in quantum mechanics, 
since any two state vectors £ Ti, such that ^ = with a S C, are in fact physically 
equivalent ~ Therefore, the appropriate quantum space of states is the space of 
rays through the origin of Ti 

v{n):=n/^, (41) 

which is called the complex projective Hilbert space" and denoted by V{T-L) for the 
infinite-dimensional T-L[38, 42]. In particular, the above-mentioned construction defines 
a vector bundle, called the "complex line bundle" , over V{7i), with C as a typical fibre, 
GL(1,C) = C — {0} = C* as the structure group and the following canonical projection 
map U -.n-^ r{n) as 

n : 1^) ^ p^, := |^)(^|. (42) 

The space V{T-L) is also called the ^quantum phase space^ which is equipped with 
the structure of a Kahler manifold {uj,g,j), induced from the corresponding structure in 
the Hilbert space, with the Riemannian structure, called the Fubini- Study metric. Its 
pull-back, according to the canonical projection map H, is defined as [36] 

:= 1 - = {d^d^) - md^)\^ (43) 

where the last equality in RHS of equation (43) holds for |<I>) = I*!*) -|- \d'^). 

In fact, the Riemannian structure of the quantum phase space, is responsible for the 
probabilistic interpretation of quantum mechanics. In other words, the geodesic distance 
between two distinguished points of the quantum phase space with respect to the Fubini- 
Study metric specifies the probability of the transition between two quantum states in 
a measurement process. As in one can show that the Schrodinger equation is the 
Hamiltonian dynamics in quantum phase space, with the additional property of being 
Killing vector field with respect to the Fubini-Study metric. 

A rather minimal approach towards the generalization of quantum mechanics is the 
notion of non-linear quantum dynamics^ associated with the extended Hamiltonian, which 
contains more general state-dependent functionals, while assuming the preservation of 
both the symplectic structure oj and the Riemmanian metric structure, i.e. the Fubini- 
Study metric g of the V{T-L). These types of generalized quantum dynamics are introduced 
by Mielnik [43], Kibble [33] and Weinberg [44], and are generally called MKW-type [45]. 

However, as it is discussed in [46], implementing the Fubini-Study metric (43) of the 
quantum phase space, does not guarantee a unique probabilistic interpretation of the 
Weinberg- type non-linear quantum mechanics [44]. The Hilbert space representation of 
the Jensen-Shannon metric (7), will be further discussed in the subsection 4.3 below. 

Accordingly, the non-linear Schrodinger time-evolution equation in % reduces to the 
non-linear Hamiltonian dynamics on V{T-L) by the projection to ViTi) of the Hamiltonian 
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vector field, which is generated by a general functional of the wavefunctions but not 
necessarily by the expectation value of the linear Hamiltonian operator. Conversely, 
the Hamiltonian flow on V{T-L) can be lifted to % which is governed by the non-linear 
Schrodinger equation. 



4.2 Hamiltonian Dynamics 

In this section, we are assuming that the non-linear quantum dynamical equation (21) is 
given. Then we show that the non-linear Schrodinger equation is in fact the Hamiltonian 
dynamics in disguise. To do so, we start from the non-linear Schrodinger equation (21) 
for one quantum particle moving in one dimension for simplicity, which is written as 

+ Vix)^ + N{p)^ (44) 



2m dx^ 



with the new Hamiltonian H = Hq + N{p) and the non-linear term as 



= Qn-Q 
(- in f 

2 



{P + Pl){p^l+ P) 



+ 



2myfp dx"^ 



(45) 



with C = h^/2mP. Since the probabilitistic interpretation of the function p{x) = 
^*{x)'^[x) is guaranteed by the continuity equation (18), the average functionals of 

denoted as Ep(//) = can easily be obtained by substituting the polar form of the 
wavefunction in 



H 



^pe (Ho + Qn{p) -Qip)) VpeT^^ dx 



AS 



which leads to 



H 



dx. 



(Mll + v + ^ln 

2m 2 



^P^{p + PiY 
{p~i + P) 



(46) 



(47) 



One has to note that the equation (47) is not an expectation value of a quantum 
linear operator, but a generalized non-linear functional of the wavefunction derived as 
the average of the extended Hamiltonian. Therefore, by using the functional derivatives 
for canonically conjugate fields p and 5, we see that the functional H is the generator of 
the Hamiltonian vector field on since 



V(p + pO" 
{p~i + p) 



-dtS 



(48) 



which is one of the Hamilton's equations of motion. The last equality in equation (48) 
comes from the real part of the polar decomposition of the non-linear Schrodinger equa- 
tion (44). For the field S one can also follow the same steps and find out that 



5H [S] {x') 
5S{x) 



— dx {pdxS) = dtp 
m 



(49) 
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which is also the Hamilton's equation of motion. As for the equation (48), the last equality 
in equation (49) comes from the imaginary part of the polar decomposition of non- 
linear Schrodinger equation. Hence, equations (48) and (49) determine the corresponding 
Hamiltonian vector field on T-L^ i.e. 

dH{.) = ^(x^,) (50) 

where r2(^, $) := 23'((^'|^>)) is the symplectic 2-form of the linear complex H in equation 
(39), and 

i! = Xfj{^):=-'jH^ = -'-{h^ + N{p))^ (51) 

with 

^^W^^ifo W + ^7V(^). (52) 

Therefore, as in the linear quantum dynamics, non-linear Schrodinger equation (44) 
can be considered as an infinite-dimensional Hamiltonian dynamics, generated by the 
generalized real- valued functional (47), in case of pure quantum states. 

Recalling the equation (42), one can obtain the projective dynamical equation and 
realize that it satisfies the following non-linear von Neumann equation 

(53) 

for pure state density operators p^, according to the x-dependence of the non-linear 
potential (45). In general, an extension from pure state non-linear Schrodinger equation 
to mixed state density operator formalism is not unique [46]. From the equation (53), it 
is implied that the scalar product for two different solutions of the non-linear Schrodinger 
equation (44), i.e. the angle between two rays corresponding to two different solutions 
of the equation (44), is not conserved. This will lead to the "mobility phenomenon", in 
the terminology used by Mielnik in [47]. 
Moreover, the following identity 

^^A^^Xjv(A^) = AX7v('f) AgC (54) 

shows that the Hamiltonain vector field onTi — {0} is homogeneous of degree one, or the 
Ep {N{p)), defined as the following 



Ep{N{p)) = J p{x)N{p)dx, (55) 

is homogeneous of degree two, i.e. 

E,(Ar(|A|V)) =|A|'E,(iV(p)). (56) 

These are the homogeneity conditions [37, 44], which imply that the non-linear dynamics 
in the Hilbert space Ti is equal to the non-linear dynamics on the unit sphere S{T-L). 
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Therefore, Ep {N{p)) is a real-valued functional, in the extent to which it contributes to 
the Ep(i^o) = (Hq) and generate a Hamiltonian flow on Ti, and also the corresponding 
Hamiltonian dynamics, i.e. the equation (53), on the quantum phase space. 

Other well-known non-linear quantum dynamics in the literature are the so called 
"non-linear Schrodinger equation" with the non-linear term as and the one intro- 

duced by Bialynicki-Birula and Mycielski in [48] with the non-linear term as —b ln(a"/9). 
In [49], Doebner et al classified the non-linear Schrodinger equations from the gauge- 
theoretical point of view based on two assumptions, 'locality' and 'separability' of the 
non-linear terms. As it is discussed in section 3, the non-linear term introduced in this 
paper, although 'separable', is manifestly 'non-local' and hence is not included in the 
gauge-theoretic classification. 

4.3 Some Comments on the Metric Structure 

The Hilbert space representaion of the Jensen-Shannon distance (7) is studied in [15], and 
it is discussed that while the analytical expression may be difficult and even impossible to 
obtain, for two close probability distributions and to the first non-vanishing term of order 
dx^, the Jensen-Shannon distance coincides with half of the pull-back of the Fubini-Study 
distance in the Hilbert space. 

In [50, 51], it is shown that constrained evolutions on (non-singular) submanifolds 
of the ViTi) are indeed non-linear quantum dynamics in general, while preserving the 
Hermitian inner product of the associated Hilbert space. Recalling the discussion of 
subsection 2.3, minimizing the Fisher metric in the statistical inference method, through 
the Lagrange multiplier technique, is equivalent with enforcing the dynamics to take place 
in a constrained submanifold of the space of normalized parametric statistical densities. 
Then the identification p{x) exp(iS/ h) maps the dynamics onto the unit sphere S{T-L) 
in complex Hilbert space. 

From geometrical point of view, by minimizing the Jensen-Shannon distance measure 
in the statistical inference method, in section 3, with the same embedding map as in the 
linear quantum mechanics, the constrained submanifold of the normalized densities from 
the space of statistical densities, i.e. A^, is mapped onto a {constrained) submanifold 
of the infinite-dimensional complex Hilbert space, which contains the solutions of the 
non-linear Schrodinger equation. 

In fact, changing the constraint by minimizing the Jensen-Shannon distance measure, 
while using the same statistical inference method and the same embedding map, can 
obviously alter the constrained submanifold in the complex Hilbert space. Considering 
the non-linear Schrodinger equation of subsection 4.2 from the constrained dynamical 
point of view, with of course different type of dynamical constraints as in the linear 
theory and the algebraic constraints discussed in [50, 51], one can conclude that the 
motion takes place in a constrained embedded submanifold of the complex Hilbert space. 
In this sense, more rigorous analyses of the embedded submanifold and the Riemannian 
metric induced on it are required, which are not studied in this paper. 
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5 Conclusions 



In this paper a non-linear version of the Schrodinger equation is derived, by applying the 
statistical inference method and using the Jensen-Shannon distance information. It is 
discussed that among all advantages it may have over the Kullback-Leibler information 
measure, such as smoothness and boundedness, the distance or metric property is the 
most significant one quantum mechanically. Therefore, it suggests that the axiomatic 
approach, discussed in [6], for using a more general information measure in the statis- 
tical inference method rather than the Fisher metric, to be equipped with the distance 
property, for the notion of distinguishability of the pure quantum states. 

The shared properties of our non-linear term with its predecessor introduced in [7] 
can be stated as follows: the scale invariance of non-linear term allows the states to be 
normalizable and the non-linear term in the Schrodinger equation is separable for either 
factorized or entangled initial states, in a system containing many particles. On the 
contrary, some of the different features with the Kullback-Leibler non-linear quantum 
dynamics are as follows: the repulsive quantum force acts on a Gaussian ansatz in all 
range of x as in the linear theory and in the perturbative regime the lowest order of 
approximation of the non-linearity term is of order of 0{P). 

Furthermore, it is shown that the non-linear quantum dynamical equation, studied in 
this paper, is in fact Hamiltonian in the complex Hilbert space, which is generated by the 
average functional of the modified quantum Hamiltonian operator as the Hamiltonian 
functional. In addition, the scale invariance of the non-linear term allows the pure quan- 
tum states to be normalizable and the structure of the unit sphere to be preserved. This 
property and the phase invariance of Kp(N{p)) are encoded in the homogeneity condition 
(56). Moreover, according to the properties of the non-linear term, the canonically pro- 
jective dynamics is in fact the non-linear von Neumann equation, for pure states density 
operators. 

Therefore, while the norms and symplectic structures of the complex Hilbert space 
are shown to be preserved by the information theoretic non-linearity introduced in this 
paper, the projective dynamics implies the preservation of the symplectic structure of 
the quantum phase space Vili.). Finally, by considering the statistical inference method 
from the constrained dynamical viewpoint, further investigations about the resulting 
embedded submanifold in % and the induced Riemannian metric on it are required. 
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